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We calculate the effects of velocity-dependent dark matter annihilation cross sections on the 
intensity of the extragalactic gamma-ray background. Our formalism does not assume a locally 
thermal distribution of dark matter particles in phase space, and is valid for arbitrary velocity- 
dependent annihilation. As concrete examples, we calculate the effects of p-wave annihilation (with 
the u-weighted cross section of av = a + bv 2 ) on the mean intensity of extragalactic gamma rays 
produced in cosmological dark matter halos. This velocity variation makes the shape of the energy 
spectrum harder, but this change in the shape is too small to see unless b/a > 10 6 . While we find 
no such models in the parameter space of the Minimal Supersymmetric Standard Model (MSSM), 
we show that it is possible to find b/a > 10 6 in the extension MSSM®[/(1)b-l- However, we find 
that the most dominant effect of the p-wave annihilation is the suppression of the amplitude of the 
gamma-ray background. A non-zero b at the dark matter freeze-out epoch requires a smaller value of 
a in order for the relic density constraint to be satisfied, suppressing the amplitude by a factor as low 
as 10 -6 for a thermal relic. Non-thermal relics will have weaker amplitude suppression. As another 
velocity-dependent effect, we calculate the spectrum for s-wave annihilation into fermions enhanced 
by the attractive Sommerfeld effect. Resonances associated with this effect result in significantly 
enhanced intensities, with a slightly softer energy spectrum. 



I. INTRODUCTION 

The simplest cosmology that consistently explains as- 
tronomical observations to date is the ACDM scenario 
where about 84% of the matter in the universe is the 
so-called dark matter [l|. All current evidence for dark 
matter is based on its large-scale gravitational effects; 
details of its fundamental nature remain elusive and are 
the center of a major campaign in modern experimental 
physics. Since the dark matter does not radiate, it must 
be non-baryonic and electrically neutral. Neutrinos can 
only account for a small fraction of the dark matter, ow- 
ing to their small masses @. Therefore, the current sim- 
plest explanation is that the major component of dark 
matter is a new stable, fundamental, massive particle 
relic from the big bang that interacts weakly with the 
standard model particles. 

Dark matter particles may be produced directly in col- 
lider experiments. They may also be detected directly in 
nuclear-recoil experiments. In addition, indirect detec- 
tion of dark matter particles from their decay products 
such as high-energy cosmic rays and gamma rays offers 
independent and promising probes of the nature of dark 
matter Indeed, there are claims that excesses of the 
high-energy cosmic-ray flux over the expected astrophysi- 
cal signals have been found: PAMELA reported an excess 
of positron fraction in the energy range of 60-100 GcV 
[J • ATIC has seen a bump in the electron plus positron 
spectrum at about 500 GeV [f| that has since been fol- 
lowed up and not seen by the calorimeter on the Fermi 
Gamma-ray Satellite Telescope @. The difficulty with 
interpreting these results in terms of annihilation or de- 
cay of dark matter particles is that there are possible 
astrophysical explanations for these observations. For 



example, there could be additional emission of positrons 
from pulsars that are not currently included in standard 
estimations of local positron abundances Q ■ There could 
also be positrons accelerated in supernova shocks around 
dying massive stars @. 

Another complication is the propagation of charged 
particles in the presence of interstellar magnetic fields 
and plasmas in the Milky Way. On the other hand, 
gamma rays do not suffer from this complication, and 
thus offer a relatively clean probe of possible dark mat- 
ter annihilation/decay signatures in the universe. In this 
paper, we describe improvements to predictions of extra- 
galactic gamma rays produced by dark matter annihila- 
tions in cosmological halos. 

To date, most of the calculations of the mean inten- 
sity and angular power spectrum of the extragalactic 
gamma-ray background have been done for s-wave dom- 
inated thermal relics with an example spectrum per an- 
nihilation motivated from the minimal supersymmetric 
standard model (MSSM) and the weak-scale annihilation 
cross section of av = 3 x 10 cm 3 /s !• (See for a 
study on the effect of the velocity-dependent Sommerfeld 
enhancement of the annihilation cross section in Galactic 
subhalos.) New methods are needed to understand the 
effects that relative velocities between annihilating par- 
ticles have on the observed gamma rays. For example, 
there are many cases in the MSSM where momentum ef- 
fects in the cross section must be pro perly implemented 
when calculating the relic density [llj . 

Instead of using the perturbation theory description of 
large-scale structure, it is more appropriate for this cal- 
culation to use the so-called "halo approach" [l2[ because 
annihilations will occur predominantly in the densest re- 
gions of space, such as the cores of dark matter halos. 
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Within this framework, we develop a new method for 
modeling the distribution of one-point mean relative ve- 
locities of dark matter particles. Our method is based 
on universal halo properties observed in the simulations 
of large scale structure that have been upheld in recent 
simulations of individual halos. 

This paper is organized as follows. Our formalism is 
described in Section [TTJ As a first application of this for- 
malism, we present a case study of p-wave annihilation in 
Section IIIII In Section IIV A[ we investigate these effects 
within the context of the MSSM. In Section HVBl we 
also consider an example of a model for which the s-wave 
component of annihilation is negligible. Although adding 
only a p-wave term to the cross section is sufficient to de- 
scribe the physics of many models, there are other inter- 
esting velocity-dependent scenarios such as Sommerfeld 
enhancements, Breit-Wigner resonances, and thresholds 
of new annihilation channels. It is conceivable that some 
of these effects will also change the spectrum per anni- 
hilation. In Section [Vj we show how these effects can 
in general be incorporated into our framework. In Sec- 
tion IVI1 we use our formalism to calculate the effects of 
Sommcrfeld-enhanced s-wave annihilation into fermions 
on the extragalactic gamma-ray spectrum. Our conclu- 
sions and the consequences of our findings are discussed 
in Section 17111 



los is well described by the Sheth-Tormen mass function 
4pj(M, z) [ID]. With this mass function, truncated halo 
profiles are used; hence, the rigid halo approximation: 
the density of matter not within any halo's virial radius 
i?vir is neglected. 

Following the analysis of [l4[, we use the virial mass- 
radius relation at redshift z of 

M = ^R 3 viT A vir p(z) (1) 

for the z-independent mean halo virial overdensity A V j r = 
180. Here, ~p(z) is the background matter density at red- 
shift z (not the critical density). Note that the Sheth- 
Tormen mass function has been calibrated against N- 
body simulations with the mass defined by the Fricnds-of- 
Friends (FoF) algorithm. This introduces an ambiguity 
as to what value of A V j r corresponds to the FoF mass. 
However, we shall ignore this subtlety because this is not 
the dominant uncertainty in our calculation. 

Although the halos in this model are not formally dis- 
joint (even with the rigidity assumption), we neglect ef- 
fects due to intersecting halos, assuming those regions to 
be rare in the model. This allows us to ignore the negli- 
gible 2 and 3-halo contributions to the 1-point statistics 
calculated below in Section HhI 



II. THE PHASE DISTRIBUTION OF HIGH 
DENSITY REGIONS OF DARK MATTER 

A. Mass function 

The spherical halo model that we adopt in this paper 
approximates the distribution of matter as an ensemble 
of universal, disjoint, rigid, spherical halos. This gives 
a good depiction of the densest regions of space, which 
are the overwhelmingly dominant sources of annihilation 
products. Although halos are ellipsoidal in general, their 
cores are seen in simulations to be distributed nearly 
spherically. 

The meaning of halo universality is that each halo's in- 
dividual properties can be expressed in terms of universal 
halo variables that are assumed to fully specify the prop- 
erties of the halo. In general, halo variables may include 
quantities such as halo position, mass, cosmological red- 
shift, concentration, formation redshift, internal angular 
momentum, etc. Universal halo functions then depend 
only on the halo variables and position within the halo. 

The statistical description promotes halo variables to 
random variables with distributions from which ensemble 
averages of the halo functions can be determined. Due 
to cosmological uniformity, statistical moments of halo 
properties are independent of halo position on a shell of 
fixed cosmological redshift. 

The simplest and most maturely developed realization 
of this framework has only two halo variables: mass, M , 
and redshift, z, of halos. The mean distribution of ha- 



B. Density profile 

The major breakthrough that makes the universal halo 
model viable is the discovery from large scale structure 
simulations that the spherically-averaged radial density 
of relaxed virialized halo structures follows a universal 
profile ph fl5| . To facilitate easier comparison with pre- 
vious results, we have decided in this paper to use the 
truncated NFW halo density profile [15( 

f — \S for r < Rvir, 

p h (r\p s ,r Sl R vir )= (2) 

[0 otherwise. 

Here, the halo profile is expressed in terms of the distance 
r from the halo center, and the halo variables are the scale 
density p s and the scale radius r s of the halo. 

It is worth noting, however, that recent simulations 
have had sufficient resolution and convergence to con- 
vincingly suggest that the NFW profile is too steep at 
T r s (l6j . Also, there are many new constraints on 
the concentration distribution from various simulations 
and astrophysical observations [TtJ ■ The robustness of 
the calculation of the extragalactic signal to these astro- 
physical uncertainties will need to be explored in future 
work. 

It is useful for what follows to explain how Eq. (|2|) is 
expressed in terms of our halo variables M and z. An 
expression for i? v ; r (M, z) is found from Eq. (JXJ) . In the 
context of truncated NFW profiles, the halo concentra- 
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tion is defined to be c = R v i r /r s , from which we get 



r s {M,c,z) = - 



3M 



47rA vir p(z) 



(3) 



Finally, integrating over the density profile for the halo 
mass and solving for the scale density gives 



Ps(c,z) 



A vh p(z)c 3 



ln(l + c) - j 



(4) 



Taking each halo to be at the mean concentration c = 
c(M, z) from [l8| gives us the expression for ph(r\M, z) 
used in the spherical halo model. 

While we shall adopt the mean concentration-mass re- 
lation of [Lsj throughout this paper, the next generaliza- 
tion of this model would be to determine a concentration- 



mass function 



dn 
dMdc 



(M,c,z). Although it appears that 



halo concentrations are distributed via a log-normal dis- 
tribution fl8| . the full joint concentration-mass distribu- 
tion needs to be determined. 



C. Velocity profile 

To insert a distribution of relative velocities in the 
spherical halo model, one would ideally want a univer- 
sal spherically-averaged halo profile of particle velocities. 
Such a profile was first observed in simulations by Taylor 
and Navarro, in 2001, in terms of the halo's density pro- 
file and its velocity variance profile a^ h (r\M, z). They 
found that 



Ph 



oc r 



(5) 



' vh 



provides an excellent fit to their numerical simulation 
[191 ] . Here, a is a constant. This result has recently 
been verified by the analysis of simulations such as the 
Aquarius simulations [l6| and GHALO (20) . which have 
much better spatial resolution. 

In general, the velocity dispersion at any position in 
the halo is not isotropic, but the simulations show that 
the anisotropy is small in the inner regions of halos, r < 
r s . As the annihilation signal is dominated by these inner 
regions, we shall assume isotropic velocity dispersions. 

Since, in this case, the density and velocity dispersion 
of a collisionless self-gravitating system is described by 
the radial Jeans equation 



d 

dr 



-r 2 d(p h <r* vh ) 



Gp h dr 



(0) 



Eq. ([5]) and Eq. ([6]) can be combined to solve for a fam- 
ily of density profiles and associated velocity dispersion 
profiles. Taylor and Navarro found a critical value for 
the proportionality constant in Eq. ([5]) that produced a 
density profile consistent with the universal profile seen 
in simulations. Dehnen and McLaughlin 21( generalized 



this analysis, generating the Dehnen-McLaughlin pro- 
files. The "critical profile" (the only physical solution 
without an outer truncation where ph becomes negative) 
occurs at a = j|, which is consistent with the values 
measured from the Aquarius simulation [l6j . The result- 
ing profile is similar to an NFW profile, but is less steep 
at the halo core. Nevertheless, it is still consistent with 
the Aquarius simulation [l6j . 

To use this result for any other density profile con- 
sistent with simulations, we can simply treat the other 
profile as an approximation of the Dehnen-McLaughlin 
solution and use Eq. ([5]) to solve for the associated halo 
velocity dispersion. We now explain how to do this for 
the NFW profile. This method can be applied to any 
other density profiles. 

When matching different profiles, it is common prac- 
tice to match them at the halo radius r_2 where the 
minus of the logarithmic slope of the density profile, 
7 = — dlnph/dlnr, is 2. However, because the power 
law profile with 7(7*) = 6 — 2a is a solution to the equa- 
tions and all other relevant solutions have a radius where 
7 takes this value, it is more convenient to parametrize 
the Dehnen-McLaughlin profiles at the radius ro where 



7(7-0) = 6 -2a. 



(7) 



Defining the halo variables po = Ph(fo) an d ct v q = 
Cvhiro), the family of Dehnen-McLaughlin profiles can 
then be parametrized by the dimcnsionless parameter 



_ 4nGp r% 



(8) 



which fixes the proportionality constant in Eq. ([5]). The 
critical NFW-like profile for a = || occurs for n = ^p. 

We match the density profiles at r$ by differentiating 
the NFW profile [Eq. ©] 



7 (r) 

and use Eq. ([7]) to find 

r (M,c,z) -- 



2r 



5 -2a 



r s {M,c,z). 



2a -3 

Putting this back into ph{To), we find 

(2a-3) 3 



(9) 



(10) 



(11) 



Finally, we use Eq. ([5]) and Eq. ([8]) to get the desired 
universal halo velocity variance profile as 



o* h (r\M, c,z) 



AttG 



pUc,z)r 2 {M,c,z) 



r (M, c,z) 



Ph {r\M,c,z) 



(12) 
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D. Derived relative velocity and mean annihilation 
cross-section profiles 

The dependence of the annihilation cross section on 
the center of mass energy of the annihilating particles 
can be written so that av is a function of the square 
relative velocity v 2 . For dark matter particles that have a 
p-wave component of annihilation, the velocity- weighted 
annihilation cross section is given by 

av = a + bv 2 (13) 

where a and b are constant coefficients and v is the rela- 
tive velocity between the annihilated p articles [22| . 

Note that, in the theory of relativistic partial wave 
analysis, the s-wave component of expansion contributes 
to both a and b [23|. However, it now appears to be 
common in the literature to simply refer to a as the s- 
wave component and b as the p-wave component. This is 
fine when the p-wave contribution to b is large compared 
to that of the s-wave, and this will be the meaning of 
those terms in the context of this paper. One just needs 
to keep in mind that, when the p-wave component of the 
partial wave is small, b will technically be dominated by 
the s-wave component of the partial wave. Nevertheless, 
the effect of b from the s-wave contribution is negligible 
in relic density and gamma-ray signal calculations, since 
it is velocity suppressed. 

Since the relevant particle physics depends on the 
square relative velocity between annihilating particles, 
we are interested in a spherically-averaged halo profile 
of the mean squared relative velocity v 2 (r) at each radial 
position within the halo. This depends on the underlying 
distribution / u (ui,r) of particle velocities Ui at position 
r. The probability distribution of square relative veloci- 
ties at that position can be calculated from 

/^(«i,r) = J d 3 u 1 d 3 u 2 / u (u 1 ,r)/ u (u 2 ,r)5(wf-|u 1 -u 2 | 2 ), 

(14) 

where S is the Dirac delta function, and the mean square 
relative velocity at that position is 

v 2 (r) = J &{v 2 )v 2 fAvl*)- (15) 

When the underlying velocity distribution / u is ther- 
mal with variance (r 2 h (r) at each halo position r, we 
find v 2 (r) = 6a 2 h (r). But halo simulations show signif- 
icant deviations from Maxwcll-Boltzmann distributions 
[HJ . Therefore, we will suppose that the universal phase 
distribution satisfies 

= Aa^(r) (16) 

for some constant A. The applicability of this relation 
should be verified by the current simulation data, which 
would give the value of A. In this paper, we shall leave 
this as a free parameter, and use A = 6 in those calcula- 
tions that require an explicit value. 



The mean relative-velocity-weighted cross section at 
each position is no longer a thermal average in general, 
and can even be taken as a universal halo function. The 
mean weighted cross section at position r of a halo is 

M„(r) = J d(«?) H(«?) r) (17) 

which, in our particular case study of Eq. (1131) . is given 
by 

[av] h (r) = a + \b(T 2 vh (r). (18) 

E. Further comments on modeling the dark matter 
distribution 

For the sample calculations in this paper, we have as- 
sumed ACDM cosmological parameters from WMAP5 [l[ 
neglecting neutrino effects: Oa = 0.721, S7b = 0.0462, 
Q c = 1 - Q A - h = 0.701, cr 8 = 0.817, and n s = 0.96. 
When calculating the halo mass function, we used the 
linear power spectrum proposed by Eisenstein and Hu 
[25| . and a critical overdensity for spherical collapse of 
1.686. 

If the Milky Way Galaxy is in an average halo of 
mass 2.0 x 10 12 M and the Solar System is 8.0 kpc 
from the center, then this model gives the local dark 
matter particle density to be 4.1 x 10~ 25 g/cm 3 = 
0.23 GeV/cm 3 , and the rms relative particle velocity 
to be 0.85 x 10~ 3 -\/A/6 c. For comparison, the lat- 
est determination of the local dark matter density 
in our solar neighborhood, based on data from both 
dynamical observations and galactic simulations, gave 
0.47 ± 0.03 ± 0.08 GeV/cm 3 ( thc nrst uncertainty is 
statistical while the second one is systematic) . This value 
takes into account the systematic effect of the dark disk 
structure which the spherical model neglects, so some un- 
derestimation by our model is expected. Thus, our model 
is able to reproduce the observed value reasonably well. 

We have made a number of additional simplifications 
while calculating thc example spectra in this paper: 

1. Thc results in our model are sensitive to the mini- 
mum mass of halos M m j n , since the cores of NFW 
halos become more dense as the halo mass de- 
creases. This low mass halo scale is determined 
by the Jeans mass, which can be determined for 
any given particle physics model [27| . These works 
show that 1O~ 12 M < M min < 10" 3 Af Q in min- 
imal supcrsymmctric extensions of the standard 
model. For the calculations in this paper, we set 
M min = 1O 6 M . Calculations using a smaller min- 
imum mass increased the intensity by a factor less 
than 10 (see Appendix |A"| . 

2. This initial work also neglects the contributions of 
halo substructures. Recent simulations have re- 
solved thc abundance of halo substructures down to 
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10 5 M Q (28|. These results are easily implemented 
in our models with recent new analytic methods for 
calculating the boost in intensity of annihilation 
gamma rays [2{|. This work is in progress. The 
halo substructure also suffers uncertainties due to 
the magnitude of Af m i n , since it also sets the mini- 
mum mass scale of sub-halos. 

3. The peculiar motions of halos, and the velocity of 
the annihilation center of mass relative to the halo 
have been neglected in this work. The correspond- 
ing corrections to our spectra due to these redshift- 
ing effects will be small. 

4. While we have included the opacity of the uni- 
verse to gamma rays in our expressions below, we 
will neglect the effect of opacity in the calcula- 
tions presented in this paper. Our ignoring the 
opacity does not change any of the conclusions of 
this paper: the suppression of photon intensity is 
quite small for photons below 10 GeV; therefore, 
its effects will only be seen for large dark mat- 
ter masses. Using the fitting formula of Steckel, 
ct al. [30[ for the opacity, we checked its effect 
for a sample model having neutralino dark mat- 
ter with a mass of 550 GeV, which is the largest 
dark matter mass appearing in the sample calcu- 
lations of this paper. The extragalactic annihila- 
tion photon intensity spectrum peaked at about 
20 GeV and the opacity began reducing the in- 
tensity more than 20% only at photon energies 
above 40 GeV, and 50% reduction at energies above 
200 GeV. For models with smaller dark matter 
masses, for instance at 150 GeV, the opacity effect 
on the gamma-ray spectrum is especially small. 



III. EXTRAGALACTIC GAMMA-RAY 
BACKGROUND: P-WAVE ANNIHILATION 

A. Formalism 

The specific intensity of gamma rays of energy E 7 from 
p-wave annihilation of dark matter particles in the direc- 
tion of n is given by @ 



7 7 (n, Ery 



dr [S p (h, r) + l] 2 [aV }^ r) W{{1 + z)E y , z) 



H{z) p[n ' Z) 



1 + - v 2 (n,z) 
a 



W((l + z)E 7 ,z) 
(19) 



where 



Mr) = & 



i 



P( r ) 



(20) 



is the dark matter overdensity. In all regions that con- 
tribute non- negligibly to the intensity, it is true that p(r), 



the dark matter density at position r, is much larger 
than p(z), the background density at rcdshift z associ- 
ated with the radial distance |r| from the observer. Here, 
[av] (r) = a+bv 2 (r) is the mean relative- velocity- weighted 
dark matter annihilation cross section at position r [3l| ; 
[<rv]o is a reference cross section that, in this context, is 
taken to be what we are referring to as the s-wave com- 
ponent of the cross section, i.e., [crv]o = lim^^o crv = a: 
v 2 (r) is the mean relative squared velocity of the dark 
matter at r; and H(z) is the Hubble function. The in- 
tensity window function W is given by 



W(E 7 ,z) 



1 



8tt 



] nl M (l + z)*^^e-^) 

(21) 

where «dm is the background dark matter number den- 
sity today, is the photon spectrum per annihilation, 
and r is the opacity of the universe to gamma rays. 

The mean intensity profile is determined from joint 
statistical moments of the density and relative veloc- 
ity fields. In our quest to quantify the new velocity- 
dependent effects, we decompose the mean intensity into 
velocity- independent and velocity-dependent terms: 



(J 7 > (£ 7 ) = <J 7 > (£ 7 ) + </ T >„ (£ 7 ) (22) 



where 



(7 7 ) (£7 7 ) = 



dz 
Wz) 



(S 2 p )(z)W((l + z)E j7 z) (23) 



is the s-wave approximation of the mean intensity (valid 
when b -C a) and 



dz 



(S 2 p v 2 )(z)W((l + z)E^z), 



(24) 

which we refer to as the p-wave component. 

Note that had we not approximated 5 p + 1 ~ S p , there 
would be an additional term with (5 p v 2 ^, but this is al- 
ways negligible when compared to the term with (5 2 v 2 }. 
These averages can be calculated in the context of the 
spherical rigid halo model by taking an ensemble average 
over the positions and masses of halos that is consistent 
with the distributions seen in the cosmological simula- 
tions. 

The statistical moments are more convenient to eval- 
uate for random variables that have vanishing mean. 
To facilitate this, we define the mean-square-relative- 
velocity overdensity 



5 V 2 (n,z) = — 1, 

v z {z) 



(25) 



where v 2 (z) is the average square one-point relative ve- 
locity at redshift z, calculated in the spherical halo model 



UhdM^(M,z)v 2 h (r\M,z), (26) 
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2 4 6 8 10 0.0 0.5 1.0 1.5 2.0 

(a) (b) 



FIG. 1: Spatially-averaged one-point velocity variance, a 2 , calculated from the spherical halo model (see Section Hi Cj) 
for the cosmology described in Section III E[ At low redshift, z< 1, the variance is exponential in z, as indicated by 
the dashed curve given by a 2 (z) ~ a 2 (0)B z with B = 0.04. The inset in the right plot shows the magnitude of the 
relative error of the exponential curve. 



or in terms of the mean velocity variance at redshift z, 



and introduced a new velocity coupling 



A 

= 4nj dM^(M,z) J dr r 2 a 2 vh (r\M, z).(27) 

A plot of a 2 (z) calculated with this model is shown by the 
solid curve in Figure [TJ This curve is not sensitive to the 
uncertainty in M m i n , since small mass halos contribute 
very little to the velocity variance. Decreasing Af m ; n from 
1O 6 M to 10- 6 M Q increases of (0) by 2x 10~ 5 % and even 
of (10) by only 1%. 

The fact that we see relative velocities increasing with 
time accentuates the fact that the particle motions are 
virial rather than thermal — increasing velocity variance 
is an indication of the growth of halo structures. In 
the model we used, the rms relative velocity today (due 
to matter in high density regions) is calculated to be 

4 



v 2 {0) = 1.4 x 10 _4 v /A/6. At low redshift, we observe 
that the velocity variance is approximately exponential 
in z, given by a%(z) ~ a%(0)B z for B w 0.04. This is 
indicated by the dashed curve in Figure Q] and shown to 
be accurate within 30% up to z = 1. Note that all of the 
calculations in this paper use the full model for cr 2 (z), 
and not the exponential approximation. 

In terms of only overdensity random variables, the 
mean intensity is written as 

xW v ((l+z)ELy,z) (28) 
where we generalized the intensity window function 

W v {E 1 ,z)=(l + -\^ v {z)]w{E 1 ,z), (29) 



1 a 



-i -l 



(30) 



that vanishes as b — > and is 1 as a — >• 0. 

After ignoring effects due to regions of intersecting ha- 
los, the statistical moments we need are simply given by 



(s 2 P ) (*) = 



A(r\M,z) 



(31) 



J dM v ' 



p 2 {z)al(z) 



(32) 



where the mass integration is taken from our minimum 
halo mass M m ; n = 1O 6 M to the maximum halo mass 
allowed by the halo concentration distribution [32[ (de- 
termined to be 2 x lO 15 A/ for the calculations in this 
paper). 

This mean intensity calculation is easily generalized for 
any relative velocity dependence of the particle physics 
model, although current derivations of the formalism re- 
quire the velocity dependence to be at least piecewise 
analytic. This will be discussed in Section IVl 



B. The contribution of the p-wave component to 
the mean annihilation gamma-ray intensity for 
arbitrary — 

The effect of the p-wave component of the cross section 
on the mean intensity of the gamma rays in this model 
is an energy-dependent intensity boost, relative to the 
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FIG. 2: Top: Aj vs. E 1 for a series of three mSUGRA 
models (at tan/3 = 10, Aq = 0, and /i > 0) having neu- 
tralino dark matter with properties given in Table |TJ For 
comparison, a power law oc E"- 018 is also plotted. Bot- 
tom: E^ (J 7 ) vs. Ej for those three models. 



s-wave approximation in Eq 

(A) (Ey) 

(J 7 } (Ely) 

where we define 



, given by: 



1 + A-Aj(£L 

a 



(33) 



A/(E 7 ) = 



(34) 

Note that this depends only on the photon spectrum per 
annihilation, the dark matter phase space distribution, 
and the opacity. 

Given that typical virial speeds for dark matter par- 
ticles at the present epoch are somewhere around uq ~ 
10 _3 c, one would only expect p-wave effects in the cross 
section to become important when - > 10 6 . Applying 
this reasoning to Eq. (f3"3"| . one would expect 



AAj ~ uTj 2 



10" 



A/ - 10" 



(35) 



Figure [2a] plots Aj(E 7 ) for some sample minimal 
supergravity (mSUGRA) models, also known as the 
constrained minimal supcrsymmetric standard model 
(CMSSM). In this paper, all MSSM particle spectra, 



cross sections, annihilation spectra, and relic densities 
are calculated using DarkSUSY 5.0.5 HU, interfaced with 
ISAJET 7.78 H3, and FeynHiggs 2.6.5.1 H|. The val- 
ues of A/ were determined numerically to within a pre- 
cision of 0.01%. The plot verifies the expected order of 
magnitude for A/, and we sec it is only mildly energy- 
dependent, slowly increasing monotonically. Over the 
energies where the intensity spectrum peaks, A/ approx- 
imately follows a power law relationship, of logarithmic 
slope 0.018, with the photon energy. The variations in 
minor structures between the curves are due to differ- 
ences in the particle spectrum per annihilation. 

After satisfying all experimental constraints, most of 
the allowed parameter space is in four regions: 

1. the focus point/hyperbolic region, with larger hig- 
gsino component in the lightest neutralino [361 ] ; 

2. the stau-neutralino co-annihilation region [37] ]: 

3. the Bulk region, where universal scalar and gaugino 
masses are small; and 

4. the heavy Higgs/A annihilation funnel, where the 
lightest neutralino mass is approximately twice the 
pseudoscalar Higgs mass [23|, [38[ . 

Among these four regions, we have given examples of 
the first three in the parameter space where tan/3 =10, 
Aq = 0, and /i > 0. Properties of those specific models 
are listed in Table [IJ Most of the Bulk region parameter 
space is ruled out by current experimental constraints; 
the Bulk region example we are presenting is on the edge 
of those constraints. 

Figure I2bl shows the annihilation gamma-ray spectra. 
Since Xb/a is much smaller than Aj 1 for each of these 
models (see the 6th column of Table [J, the velocity term 
of the intensity contributes negligibly to these spectra. 
However, the p-wave strength b/a of the cross section 
does have consequences for the normalization of the in- 
tensity curve, especially for the co-annihilation region. 
We shall explain this in the next subsection and in Sec- 
tion HO] 



C. P-wave suppression effect on the amplitude of 
the mean annihilation gamma-ray intensity 

The relic density of dark matter has reached the sta- 
tus of being an important constraint of particle physics 
models, owing to the fact that the magnitude of the relic 
density in our universe is precisely determined [l[, and 
that the technology needed to calculate the relic density 
for any particle physics model is mature. 

In the scenario of a thermally produced dark mat- 
ter relic, numerical solutions of the Boltzmann equation 
present an approximate picture that holds, up to nu- 
merical corrections, typically within an order of magni- 
tude: dark matter particles interact in thermal equilib- 
rium with the big bang plasma until the thermal tem- 
perature is too low to produce new dark matter particles 
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(xlO cm /sj 


Focus Point 


2569 


395 


150 


1.9 


1.8 


0.114 


22.9 


2.6 


Bulk 


79 


171 


62.3 


0.27 


57.5 


0.114 


22.5 


3.8 


Co- annihilation 


79 


373.7 


150 


0.0019 


378.8 


0.113 


24.0 


5.8 



TABLE I: Sample mSUGRA models with parameters tan/3 = 10, A n = 0, and fi > 0. The first two columns show 
the input model parameters. The other columns show the calculated dark matter properties that are relevant to our 
discussion. The last column shows the thermally-averaged av at freezeout, including co-annihilations. 



and the number density of dark matter particles is low 
enough so that annihilations are rare [H . After this time 
of freezeout, the number of dark matter particles in the 
universe per comoving volume is essentially constant (ex- 
cept for some residual annihilation occurring at very low 
annihilation rates). 

Therefore, the relic density magnitude can be ac- 
counted for, except for some important exceptions, by 
two quantities: the freezeout temperature normally ex- 
pressed as Xf = jtidm/T/, and the mean value of the dark 
matter annihilation cross section at freezeout [crz;]^. To 
thermally produce the correct relic density, the approxi- 
mate value of the velocity weighed annihilation cross sec- 
tion needs to be [uv]f w 3 x 10~ 26 cm 3 /s. 

One important exception to the above argument is the 
following: if co-annihilations are present, then the an- 
nihilation cross section at freezeout is replaced by the 
larger effective co-annihilation cross section and the re- 
sulting relic density is lower. For precision, we take the 
definition of Xf to be at the thermal temperature where 
the dark matter number density is twice its thermal equi- 
librium value. 

When in thermal equilibrium, the mean square relative 
velocity of the dark matter is related to the thermal tem- 
perature by v 2 = Q/ Xf. Then, insofar as the approximate 
picture is valid, to meet the dark matter relic density con- 
straint, a theory with p-wave annihilation strength b/a 
must have an s-wave component of approximately 



to satisfy the relic density constraint. Here, when cal- 
culating Aj in this expression, we use the cross section 
at freezeout, [cv]f, as the reference cross section [av] 
appearing in the window function given in Eq. (|21[) . 

This effect is plotted for mild values of b/ a in Figure l3"al 
for x f = 24 and AA/ = 5 x 1CT 7 . Since AA/ < 6/x f , 
there are 3 regions of interest as can be seen explicitly in 
Figure |3b] 

1. When b/a <C Xf/6 ~ 4, the theory is s-wave dom- 
inated, and the p-wave contributes very little to 
both the relic density calculation and the annihi- 
lation photon intensity spectrum. We find when 
b/a < 1, the intensity suppression is linear in b/a 
with slope —Q/xf ~ —1/4. 

2. When x//6 <C b/a <C (AA/) -1 , it is important to 
include the p-wave for the accurate calculation of 
the relic density, resulting in a suppression of the 
s-wave of the annihilation cross section. But the 
p-wave contribution to the shape of the intensity 
spectrum is negligible in this region. 

3. When b/a > (AA/) -1 ~ 10 6 , the theory is p-wave 
dominated and the s-wave component contributes 
little to both calculations. In this regime, p-wave 
suppression is maximal with a suppression factor of 
X/AA//6 ~ XfAj ~ 10~ 6 relative to an equivalent 
s-wave dominated theory. 



1 



(36) In summary, for each of the three regions, the p-wave 
suppression goes like 



in the absence of co-annihilations. That is, the pres- 
ence of p-wave annihilation requires the theory to have 
a smaller s-wave component in order to satisfy the relic 
density constraint. This is the p-wave suppression effect. 

If we were to compare the mean intensity of extra- 
galactic annihilation photons with energy E 1 in an s- 
wave model with b/a = to another model with the 
same annihilation spectrum, dark matter particle mass, 
and freezeout temperature, but with p-wave strength b/a, 
then the intensity due to a thermal relic would need to 
be p-wave suppressed by an approximate factor of 



(I, 



IS) 



(7 7 )(£y|0) 



l + (^)AA/(£ 7 ) 
1 + m^ 

V a / Xf 



(37) 



1 — — (-) £ < ^ 

x f \aJ a ^ 6 ' 




^ « I < [AA/C^)]- 1 , 

I » [XAj(E 7 )}-\ 

(38) 

If at freezeout there are co-annihilations that increase 
the effective dark matter annihilation cross section, they 
will contribute further suppression factors to the cross 
section today, and thus to the mean intensity spectrum. 

We conclude that, although the p-wave component 
does not contribute strongly to the mean intensity of the 
cxtragalactic annihilation gamma rays unless b/a > 10 6 , 
even mild values of b affect the relic density constraint, 
resulting in significant reductions to the intensity today. 
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FIG. 3: Top: The p-wave suppression factor for the mean 
intensity of extragalactic dark matter annihilation pho- 
tons [Eq. (|3"7|) ] as a function of b/a, for typical values of 
xj = 24 and AAj = 5x 1CP 7 . The suppression is shown 
here for small values of b/a. Bottom: The same plot, but 
on log scale and for a wider range of b/a. 



Of course, the relic density constraint does not apply to 
non-thermal relics, such as the post-freezeout reheating 
scenario 1391 . 



IV. EXAMPLE P-WAVE EFFECTS IN 
SUPERSYMMETRIC EXTENSIONS OF THE 
STANDARD MODEL 

A. P-wave annihilation in the MSSM 



In Section [ill Bl we justified the intuitive notion that, 
in order for p-wave annihilation effects to be important 
when calculating the mean intensity spectrum of ex- 
tragalactic dark matter annihilation photons, one needs 



b/a > 10 6 . In Section llll CI we argued that this is true, so 
long as the correct associated s-wave annihilation com- 
ponent needed to satisfy the relic density constraint is 
used. 

In the case of the MSSM, we find model parameters 
that satisfy the relic density constraint and determine 
the s-wave and p-wave components of the dark matter 
relative- velocity- weighted annihilation cross section for 
that theory [23]. It is interesting to consider which re- 
gions of the MSSM have the largest ratio of p-wave to 
s-wave components. 

The first step is recognizing that ncutralino dark mat- 
ter has certain annihilation channels where the s-wave 
components are helicity-suppressed. For annihilation 
into a fermion anti-fermion pair, the s-wave annihilation 
component is suppressed by (m//?n^o) 2 , the squared ra- 
tio of the fermion mass to the neutralino mass. This in- 
cludes t and u-channel sfcrmion exchange, and s-channcl 
mediation by the Z boson or neutral Higgs bosons. Each 
contribution to the total cross-section amplitude has an 
approximate factor of (mjv[/m^o)~ 2 , where rriM is the 
mediator mass. This factor enhances the channels me- 
diated by the Z and lightest Higgs, but suppresses the 
sfcrmion and heavy Higgs channels. 

In the parameter space where the neutralino is very 
nearly pure bino and the magnitude of the Higgs super- 
potential coupling n is much larger than the first soft 
gaugino mass Mi, the annihilation is dominated by the 
sfcrmion exchange. This allows us to greatly reduce the 
contributions of Z or Higgs mediation, as well as all other 
annihilation channels that are not helicity-suppressed. 

In the case that the sfermion masses are unified at the 
GUT scale, the heavy third generation fermion channels 
dominate the s-wave annihilations, due to their larger 
masses compared with the other fcrmions. In parameter 
space where the neutralino is nearly pure bino, squark 
masses are much larger than sfermion masses. However, 
the b quark is sufficiently massive when compared to the 
t lepton that it still has a significant branching ratio. 
Also, remaining modest amounts of Z mediation will add 
to the bb production. The production of t quarks be- 
comes important if the neutralino is massive enough to 
kinematically allow it. 

We can increase b/a further by taking advantage of 
the sfermion mass suppression of the cross section and 
considering large non-universal 3rd generation sfermion 
masses at the GUT scale. In this parameter space of the 
MSSM, where the ncutralino is nearly pure bino and an- 
nihilation into 3rd generation fcrmions is suppressed, we 
would expect the dominant contribution to the s-wave 
component of the cross section to be proportional to 
(m^/m^o) 2 < 1CP 6 , relative to the p-wave component. 

However, there are loop processes that become dom- 
inant at this point, generated when the fermion anti- 
fermion pair close the loop and two gauge bosons come 
off the internal lines [4(|, such as two photons, photon 
and Z, or two gluons. The amplitude due to these loops 
provides a hard lower bound on the s-wave component 
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of the cross section on the order of 10~ 29 cm 3 /s, keeping 
b/a < 10 4 in the MSSM parameter space that satisfies 
the relic density constraint. Correspondingly, the p-wave 
intensity term will have a magnitude of less than a per- 
cent of the s-wave approximation. But as we explained 
in Section IIII CI such large values of b/a require a sig- 
nificant reduction in the s-wave component of the cross 
section in order to satisfy the relic density constraint. 

Therefore we conclude that, in the MSSM, the s-wave 
approximation of the intensity calculation in Eq. (|23j) 
gives very accurate results in the MSSM, as long as the 
correct annihilation cross section is used. We find that 
p-wave suppression of the s-wave component is very sig- 
nificant in parts of the MSSM parameter space. The 
intensity is further suppressed in parameter space with 
significant co-annihilations at dark matter freezeout. 

As an example, we showed in Figure [2] a comparison 
of the predicted extragalactic annihilation spectrum for 
three different MSSM models in the mSUGRA parame- 
ter space. We expect conclusions about all other regions 
to follow similarly from the results of this section, under- 
stood from Figure [3bl Although we see slight variations 
in the models' spectral shape due to differences in the 
photon spectrum per annihilation in the different regions, 
the dominant effect is the normalization of extragalac- 
tic intensity. The co-annihilation region spectrum is so 
strongly p-wave suppressed (and also mildly suppressed 
further because of the co-annihilation factors) that it is 
a factor of 0(1O -2 ) less intense than the corresponding 
focus point/hyperbolic region spectrum that has nearly 
no p-wave or co-annihilation suppression. 

The large p-wave strength of b/a = 378.8 of the co- 
annihilation region can be accounted for by the fact 
that the lightest neutralino is about 99% bino and 1% 
higgsino. Therefore, its annihilation is dominated by 
helicity-suppressed processes: 66% into bb and 24% into 
t + t~. In contrast, the focus point/hyperbolic region's 
lightest neutralino is 72% bino, 26% higgsino, and 2% 
wino. It annihilates predominantly into W + W~ (54%), 
bb (20%), and ZZ (18%). 

Note that not all co-annihilation regions of mSUGRA 
are strongly p-wave suppressed. At larger tan/3, the s- 
wave annihilation into fcrmions mediated by an s-channel 
pscudoscalar Higgs is enhanced, producing mostly bb 
pairs, and the p-wave strength decreases back down to 
0(1) in the co-annihilation region. For example, the up- 
per curve in Figure 2] shows the intensity spectrum for 
a co-annihilation region mSUGRA model as before, ex- 
cept tan /3 = 50 and m^a = 550 GcV. In this model, 
the p-wave strength of b/a = 4.8 is lower because of the 
contribution of the pseudoscalar-mediated annihilation 
channel, even though the lightest neutralino is 99.8% 
bino. Near the edge of the current b — > 57 bounds 
on the co-annihilation region parameter space, where 
m^o = 183 GeV, the p-wave strength is even lower at 
b/a = 2.7. The p-wave strength increases for higher 
neutralino mass in this co-annihilation region because: 
\x increases rapidly with neutralino mass, and the pseu- 
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FIG. 4: Top: The lower curve is the intensity spectrum 
for a pure p-wave process. This scenario is realized in 
the MSSM(8)(7(1)_b_l with right handed sneutrino dark 
matter. The sneutrino, here of mass 550 GeV, annihi- 
lates through a Z' resonance into fermion anti-fermion 
pairs. The Z' for this plot had mass 1300 GeV and width 
17 GeV. For contrast, the upper curve is due to annihi- 
lations of a 550 GeV neutralino in the co-annihilation 
region of mSUGRA with tan/3 = 50, A§ = 0, and 
fi > 0. This model has a relatively mild p-wave strength 
of b/a = 4.8. Bottom: The associated A/ for the two 
models, given by Eq. (|34l) . For reference, the same power 
law as in Figure [5a] is also shown. 



doscalar mass increases with /i; therefore, the s-wave an- 
nihilation with pscudoscalar mediation becomes impor- 
tant. 



B. A p-wave dominated scenario 

Can we find an extension of the MSSM which would 
yield a larger value of b/a? 

One interesting extension of the MSSM is that of 
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adding a U(1)b-l g au g c symmetry, based on a charge 



of baryon number minus lepton number [41[ . This model 



is interesting because it provides a gauge structure that 
accommodates the right-handed neutrino. 



V. INTENSITY SPECTRA OF 
EXTRAGALACTIC ANNIHILATION GAMMA 
RAYS FOR GENERAL VELOCITY-DEPENDENT 
PARTICLE MODELS 



In this model, the lightest supersymmetric particle 
(LSP) can be the lightest right sneutrino or the lightest 
neutralino. Here, we consider a parameter space where 
the right sneutrino is the LSP. Thus it is neutral to stan- 
dard model charges and the only gauge interaction is with 
the Z' (and B — L neutralinos) via its lepton charge. It 
also interacts with the B — L Higgs fields via the £>-term. 
Possible s-wave annihilation processes are annihilation 
into neutrinos mediated by B — L neutralinos, and into 
B — L Higgs via s-channel Z' or Higgs fields. However, 
there is parameter space where the B — L Higgs' and 
neutralinos are massive compared to twice the sneutrino 
mass, and can be neglected in this discussion. In this re- 
gion, sneutrino annihilation is (at tree level) exclusively 
s-channel via the Z' into fermion anti-fermion pairs. In 
this process, s-wave annihilation is completely forbidden. 
So, this is an example of a pure p-wave annihilation pro- 
cess. 

At one loop, an s-wave component is generated, but is 
strongly coupling-suppressed when compared to the tree 
level p-wave cross section, with additional factors of a 4 
or g' 4 a 2 , where a is the fine structure constant (or the 
strong force constant in the case where gluons are emitted 
from quarks, instead of photons) and g' is the U(1)b-l 
gauge coupling, which we take to be 0.4. Thus, we would 
estimate b/a in this scenario to be > 10 8 , completely 
p-wave dominated. We carry out a calculation of the 
intensity spectrum due to annihilations at the tree level, 
neglecting the small s-wave component altogether. 

For our example, we consider a model with spec- 
trum m = 550 GeV, M z > = 1300 GcV, and Z' width 
T z' = 17 GcV. Here, the sneutrino is at a near resonance 
with the Z' , which allows the relic density constraint to 
be satisfied. The photon spectrum per annihilation was 
simulated with Pythia 8.135 [42|. The resulting intensity 
spectrum is shown in the bottom curve of Figure 2] The 
low intensity in this model shows the expected full p-wave 
suppression of 10 6 , as described in Section Till Gl Here, 
the differences in spectrum shape between the two mod- 
els will be due to both differences in photon spectrum per 
annihilation, and the fact that the lower intensity curve 
is directly coupled to the current dark matter velocity 
distribution, whereas the upper curve is not. 

We note that in the context of the B — L extension 
of the MSSM, there are also parameter regions with sig- 
nificant s-wave annihilation components. For example, 
regions with light B—L Higgs fields open up sneutrino an- 
nihilation channels involving only scalar fields whose con- 
tribution to the cross section is momentum-independent 
0. 



Adding a p-wave term to the annihilation cross section 
Eq. (fT3"| is sufficient to describe the physics at energies 
of the current dark matter distribution for many particle 
physics models. However, there is a variety of additional 
momentum-dependencies that can arise. 

It is therefore useful to develop a formalism for calcu- 
lating mean extragalactic annihilation photon intensity 
spectra for general velocity-dependent particle physics 
models. Given a relative-vclocity-wcighted annihilation 
cross section as a function of squared relative velocity 
[cri;](i; 2 ), the mean photon intensity due to extragalactic 
annihilations is simply 



(7 7 ) (£ 7 ) = 



dz 



^^W((l + *)£,,*), 

p [Ho / 

(39) 

where, again, [crv]o is any convenient reference cross sec- 
tion used in the intensity window function W, in Eq. (j2"T|) . 

To calculate the ensemble average, we use the fol- 
lowing result from the spherical halo model. Let 
Xi(r, z), . . . , Xjv(r, z) be any fields associated with uni- 
versal halo functions X\h{r\M, z), . . . , X]yh(r\M, z). The 
overdensity of Xi is 



where 



Xi(r,z) 



dn 



(40) 



Xi(z) = j dM—(M,z) J d^X, h (r\M,z). (41) 

The ensemble average of any product of these overdensi- 
ties, all evaluated at one position, at redshift z is 



n 71 

dM y ' ' 



N 



X lh (r\M, z) 



i=l 



Xi(z) 



(42) 

where we have used the disjoint, rigid halo approximation 
[44]] . Using Taylor's theorem, we find a formula for the 
ensemble average for any piecewise analytic function F of 
the overdensities of quantities associated with universal 
halo functions. 

(f(5 Xi ,...,S Xn ))(z) 

]3 ^fX lh (r\M,z) X Nh (r\M,z) 



Xi{z) 



X N {z) 



(43) 



In the case of Eq. (|3T))) . we re-express av as a function 
of 5 V 2 at redshift z by substituting v 2 = Xa 2 (z)(l + 5 V 2) . 
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We then obtain 



r /<!,,(r|M,z) 



(44) 



There are calculations where it may be more convenient 
just to, when possible, treat the cross section itself as a 
universal halo function, given by Eq. (|17|) . In those situa- 
tions, one can define the mean cross section av at rcdshift 
z and associated overdensity S av field in the usual way, 
with Eq. (|4"0j) and Eq. (|4"Tj) . Then we have the alternative 
expression 

av(z) 



with 7710 the mediator mass and mx the dark matter 
mass. For convenience, we also define e v = ^. 

Following the derivation by [46j], 5 is extracted from 
the solution of the Schrodinger equation 

$" + 1$' + (1 - {/)$ = 

with boundary conditions $(0) = 1 and $'(0) = —l/e v , 
and where 

is the normalized potential. This is more illuminating 
with ip( x ) = 2G*l>(a:), in which case the Schrodinger equa- 
tion becomes 



CT?7(z) 



3 p 2 Jr|M,z) 
P 2 W 



(45) 



Similarly, if the annihilation product spectrum per an- 
nihilation should vary significantly over the range 
of annihilation center-of-momentum energies that occur, 
then it may also be expressed as a function of relative 
velocity and included in the ensemble average. 



VI. SOMMERFELD ENHANCEMENT AND 
RESONANCE EFFECTS IN THE MEAN 
INTENSITY SPECTRUM 



The Sommcrfeld enhancement [4 51 - 44 7j of dark matter 
annihilation occurs in cases where the annihilation is me- 
diated by an attractive Yukawa force through a scalar 
or vector boson. Non-perturbative resonant boson ex- 
change between annihilating particles is demonstrated to 
result in a significant enhancement of the cross-section 
that grows as relative particle motion decreases. 

Although this enhancement occurs for any partial wave 
plrl Wi\ , we will (for simplicity) present intensity spectra 
for s-wave Sommcrfeld enhancement. It is expressed in 
the form 



[<rv](v) 



S (l\ e *) 



av\ 



(46) 



where [crv]o is the relative- velocity- weighted annihilation 
cross section at tree level, assumed constant at all rele- 
vant energies for this model, and set for our calculations 
to [av]f = 3 x 10 26 cm 3 /s to satisfy the relic density (al- 
though, in careful calculations, the enhancement effect 
can have some affect on the relic density calculation [48| ) . 
Here, a is the Yukawa coupling between the dark matter 
and mediator, and 



4" + (l - U)tp = o 



(47) 



with the necessary solution near the boundary of 
lim^-^o "0 — x ~ x 2 / e v and lim 2; _ s .o ip 1 = 1 — 2x/e v . 

It is now easy to see that tp converges very quickly to 
a sinusoid as x increases. S is simply the inverse square 
of the amplitude of ip f ar from the origin. One can in- 
tegrate Eq. (|4T|) to x = xm large enough that U{xm) is 
sufficiently negligible and simply evaluate 



S = 



1 



i/j 2 (x m ) +V" 2 (x m ) 



(48) 



Figure [5a] shows the enhancement for e v = 1, 0.1, and 
0.01. As relative velocity decreases, a series of Sommer- 
feld resonances reveals itself. Let us define the locations 
of the resonances to be at = for n = 1, 2, Ana- 
lytic approximations show the first few s-wave resonances 
to be near e? ~ 6/(n7r) 2 . Also, 



A 



lira S{e v \e2 ) 



where A i ' are constants. We find the first resonance 
at = 0.595 306 210 530 309 to have A W = 3.37286 
up to the given precisions. Away from the resonances, 
the enhancement saturates to a constant value as v di- 
minishes. Until the resonances arise at low e^, S follows 
the Coulomb-Sommerfeld enhancement in the center-of- 
momcntum frame 



S{e v \Q) 



2n/e v 



1 



-27r/e„ 



mx ol 



which goes like e v 1 at low v. Examples of the relative 
velocity dependence of each of these cases are shown in 
Figure [5b] for a = 0.01. 

Figure|H]shows the results for the extragalactic gamma- 
ray intensity due to annihilations for theories with these 
cross sections. Again for simplicity, we consider theo- 
ries that annihilate exclusively into lepton anti-lcpton 
pairs. The photon radiation spectra per annihilation 
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FIG. 5: Left: The s-wave Sommerfeld enhancement S(e v \e ( j > ) vs. for e v = 1, 0.1, and 0.01 from bottom to top, 
respectively. The vertical lines correspond, from right to left, to = 1.107, 0.635, and the first Sommerfeld resonance 
e^ 1 ' described in the text. Right: S^le^) vs. v 2 for a = 0.01. The solid curves show the enhancement for the same 
three values of specified earlier from bottom to top, and the dashed curve shows the Coulomb case where = 0. 




FIG. 6: Extragalactic annihilation gamma-ray intensities for s-wave theories, sample Sommerfeld-enhanced s-wave 
theories, and Sommerfeld s-wave resonances. For these models, we set a = 0.01 and mx = 200 GeV, and used A = 6 
when generating [crv](<5„2) in Eq. (j44[) . The bottom dotted curves show the intensity for no Sommerfeld enhancement, 
the solid lines show the Coulomb-Sommerfeld resonance, the top dot-dashed curves show the first Sommerfeld reso- 
nance, and the long (short) dashed curves show the Sommerfeld enhanced intensity where the enhancement saturates 
at 10 (1000). Left: primary photon radiation from annihilation into electron-positron pairs. Right: annihilation into 

T + T~. 



for these processes were also simulated with Pythia. In 
these examples, the dark matter mass is taken to be 
mx = 200 GeV. Higher masses were explored; they sim- 
ply gave the same results, scaled up in energy. Annihi- 
lation into /i + fi~~ is visually indistinguishable from the 
electron case, with a slightly higher intensity. The ratio 
of each enhanced case to its respective s-wave approxi- 
mation (where S = 1) is shown in Figure [7] Here we find 
that, for a = 0.01, the intensity enhancement is nearly 
uniform over most of the Sommerfeld parameter space. 
This follows from the fact that, as we see in Figure I5b[ 
the cross section has already saturated at the relative 
velocities important today, around 10 -4 , unless you are 
very, very close to a resonance value of e^. However, for 



smaller values of a, the graph in Figure l5bl shifts to the 
left and the cross section may not necessarily be com- 
pletely saturated today for larger enhancements. This 
results in a smaller cross section at high energies and has 
the effect of widening the intensity peak and shifting the 
maximum to smaller energy This is what is observed 
with the resonance cases: = (an example of a v~ 1 
resonance), and = ' (i>~ 2 resonances). 

To see the detail of the variation of (7 7 ) / (/ 7 ) at peak 
intensities, plots of each ratio on a blown up linear scale 
are shown in Figure [8] For comparison, the enhance- 
ment at saturation for = 1.107 was 10.00, and for 
£0 = 0.635 it was 1004. The intensity ratio for annihila- 
tion into taus in the saturated examples was very similar 
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FIG. 7: Ratios of the predicted Sommerfeld enhanced intensities to the unenhanccd intensity. Left: primary photon 
production from annihilation exclusively into electron-positron pairs. Right: annihilation into t + t~. 
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FIG. 8: The intensity ratios in Figure [7] The solid curves are for annihilation into e + e , and the dot-dashed curves 
show the results for r pair production. 



to the electron-production results. However, we can see 
that differences in the spectrum per annihilation become 
important at the resonances. 

In principle, the resonance cross-sections can break 
unitarity bounds for s-wave annihilation [49| . If the scat- 
tering operator conserves angular momentum and is uni- 
tary, then the weighted s-wave annihilation cross section 
must satisfy [5(| 



av < 



4tt 

m 2 x v 



(49) 



In the low-w limit v <C a for the Coulomb case, this pro- 
vides an upper bound on the mass for a given coupling: 



(30TeV) 2 , fore = O (50) 



for the value of [cv]f in our model. For the first Som- 



merfeld resonance, we require 



> 



A 



4tt 



[av] f 



mxOi 
40TeV 



for e„ 



(51) 

As long as the enhancement saturates before getting to 
relative velocities below v m - m , then the theory is consis- 
tent. Otherwise, some neglected mo del- dependent effects 
(such as finite widths or non-pcrturbative dynamics) be- 
come important and must be taken into account. A 
model near resonance that saturates below scales that 
contribute to the intensity would be indistinguishable 
from the resonance intensity, but would still satisfy uni- 
tarity constraints. 

We finish this section by mentioning another common 
example of resonance effect: the Breit-Wigner resonance 
due, for example, to the mass of the mediator being at the 
energy of the annihilating particles. At non-rclativistic 
speeds, this phenomenon also results in saturated cross 
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sections, but in certain situations such as when the dark 
matter mass is very near half the resonance mass and 
much larger than the resonance width, we have av ~ v~ A 
for non-relativistic relative velocities above the satura- 
tion scale (see Appendix [Bj . 



VII. DISCUSSION 

In this paper, we have studied how velocity-dependent 
annihilation of dark matter particles can affect the 
mean intensity spectrum of primarily produced, extra- 
galactic, annihilation gamma-ray background. Velocity- 
dependence of dark matter annihilation contains infor- 
mation about the spin of the dark matter and its anni- 
hilation mediators, and about various possible resonance 
effects. 

We first considered the effects of p-wave annihilation, 
where the annihilation cross section is well approximated 
by av = a + bv 2 . This is an important feature in min- 
imal supersymmetric extensions of the standard model 
where the s-wave component can be strongly helicity- 
supprcssed. If the relative- velocity- weighted cross section 
is quadratic in relative velocity at all energy scales up to 
freezcout, then the relic density constraint requires the 
s-wave component of the cross section to be reduced for 
larger p-wave components, resulting in a suppression of 
the extragalactic gamma-ray intensity today. 

Other than setting the scale of the s-wave annihilation 
component, we find that the p-wave does not contribute 
to the intensity unless the coefficient b is at least 6 or- 
ders of magnitude greater than the s-wave component. 
We do not find such a scenario within the MSSM (we 
found MSSM p-waves only as high as 4 orders of mag- 
nitude greater than the s-wave component); therefore, 
differences in mean intensities in the MSSM are due to 
factors from different s-wave components (from the pres- 
ence of p-wave annihilation or co-annihilations), or major 
differences in the photon spectrum per annihilation. 

Although there was no model found in the MSSM 
where p-wave annihilation dominates today (with b/a > 
10 6 ), we presented such a scenario within the context 
of MSSM<E){/(1)b_l, where baryon number minus lepton 
number is an additional gauge charge with new Z' vector 
boson. This model contains sneutrino dark matter at a 
near resonance with the Z' in order to satisfy the relic 
density constraint, but the cross section remains propor- 
tional to the squared relative velocity at tree level. In 
this case, the intensity spectrum is maximally coupled 
to the velocity distribution, but the dominant feature is 
again the strong 10 -6 suppression in the normalization 
of the intensity curve. 

However, one cannot conclude that theories with 
strong p-wave annihilation will always result in an un- 
observably small extragalactic gamma-ray signal, since 
non-thermal freczeout of the dark matter can allow for 
larger intensities. 

A possible cross section feature that we found to en- 



hance the mean intensity of extragalactic annihilation 
gamma rays was the Sommerfeld-enhanced s-wave cross 
section. If the dark matter attracts itself, then s-wave 
annihilation becomes stronger at low energies. Again, 
the dominant feature is a change in the normalization of 
the mean intensity (J 7 ) (E^), but this time an increase by 
as high as 10 4 . The Sommcrfcld enhancement also has 
the effect of shifting the peak of E 2 (J 7 ) to lower pho- 
ton energy E 1 , thereby broadening the intensity peak. 
The Sommerfeld resonances provide an example of v~ 2 
behavior for the velocity- weighted cross section, or 
behavior for an extremely light or massless mediator. An 
example of i>~ 4 behavior can be found in Breit-Wigner 
resonances for certain ranges of v. 

In all of the cases for which wc calculated intensities 
and compared them to an associated constant av model, 
we find that the dominant effect is a change in the nor- 
malization of the mean intensity. By considering veloc- 
ity effects alone, we could account for modification in 
intensity magnitude over 10 orders of magnitude (up to 
6 orders decrease from p-wave effects, and over 4 orders 
increase possible from Sommerfeld enhancements). Non- 
thcrmally produced dark matter relics complicate this 
even further. 

It is straightforward to extend our formalism to the cal- 
culations of the anisotropy in the gamma-ray background 
Different particle physics models will have unique 
angular power spectral features, independent of inten- 
sity normalization [Hl|. For example, in a Sommerfeld- 
enhanced theory, small halos have smaller virial motions 
and will appear slightly brighter than for a theory with- 
out enhancement. 

Future work will also need to explore the robustness 
of these calculations to the astrophysical uncertainties, 
including the universal halo functions, and distribution 
of halos and subhalos. 

If the LHC is to detect a candidate dark matter parti- 
cle, the discovery will require verification of consistency 
with the existing dark matter population through direct 
and/or indirect observations. The methods outlined here 
are able to improve predictions of indirect signals for cer- 
tain particle models and expand the reach of interpreta- 
tions of analysis of the observed signals. In order to be 
able to extract new particle physics and astrophysics in- 
formation from the intensity of an indirect signal, the 
challenge remains to identify robust features that can be 
disentangled from the uncertainties of the dark matter 
intrinsic properties and astrophysical distribution. 
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FIG. 9: Intensity spectra of dark matter annihilation 
gamma-rays for the MSSM focus point model described 
m Section HITbI with minimum halo mass M m i n = 
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the minimum halo mass from 10 b A/q to 
0, for the focus point and bulk particle mod- 
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els described in Section MI Bl The bulk model 
is the curve over slightly lower energies, since 
it has a less massive dark matter particle mass. 
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Appendix A: Effect of Low Mass Halos 

The intensity of photons from extragalactic annihilat- 
ing dark matter is sensitive to the scale of minimal halo 
mass, M m i n , which is dependent on details of the dark 
matter self-interaction. For simplicity, the spectra cal- 
culated in this paper set Af min = 1O 6 M . To gam un- 
derstanding of the effect of low mass halos, we plot in 
Figure [9J the intensity for a model with different values of 
Afmin- The intensity enhancement from M m - ln = 10 6 A/q, 
at the spectral peak, is 4.6 for M min = 1O~ 6 M and 7.3 
for M min = 1O" 12 M . 

Figure [TU] shows the enhancement factor for M m - m = 
1O _6 M (the curve offset to the right). We see the en- 
hancement is essentially uniform, but is slightly weaker 
at the spectral peak than at the tails. For comparison, 
Figure [TU] also shows the plot for the bulk model, also 



described in Section IIII Bl We see that decreasing the 
minimum halo mass has practically the same effect for 
the different models. 



Appendix B: Breit-Wigner Resonances 

If the dark matter annihilates into an unstable particle 
of mass M and decay rate T <C M that then decays, and 
the center-of-momentum energy E of the annihilation is 
near M, then the annihilation cross section is of the form 

m 



If we define 



1 



M 2 T 2 



E^JE 2 - Am 2 x (E 2 - M 2 ) 2 + M 2 T 2 ' 



(Bl) 



M 
M 

Am — 

2 

Am 
~W 



mx, and 



Am 



1 



Am 



then, for any relative velocity v of the annihilating parti- 
cles and spectral separation Am, the velocity dependence 
of the cross section near the resonance is 



,(v) oc 



i 



-4Am 



r i 



-(f) 2 



(B2) 



For a non-rclativistic resonance, v <C 1 and |Am| <C M 
giving 



(v) 



1 



4Am 



(B3) 



where, in this context, Am = Am/M and [av] r is the 
velocity-weighted cross section at the resonance energy 
where v = v r = 4V Am (which is unattainable in the 
case that M < 2m). 

If there are no other significant additional features in 
the cross section from the resonance scale v r to the freeze- 
out scale Vf = 6Tf/mx ~ 1/2 for freezeout temperature 
Tf, then [av] r can be correlated to the mean cross section 
at freezeout [crv]f ~3x 10~ 26 cm 3 /s: 

r M/ 



(i6r) 2 

where we assumed Vf ~ 1/2, |Am| -C 1/64, and T <C 
1/16, and we neglected some 0(1) constants due to Vf 
being near 1. In this kind of theory, 

n 5/2 



av 



,(v) 



(16) 2 



1 - 



4Am 



r 2 



(B4) 
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up to the freezeout scale where r\ is an 0(1) constant 
that takes into account our approximations of the relic 
density calculation and relativistic freezeout velocities. 

The velocity dependence for s-wave annihilation via 
a non-relativistic (v r <C 1 and |Am| -C M/64) Breit- 
Wigner resonance of small width (r <C M/16) is therefore 
found to have the broad behavior of 



\<JV\ 



,(«) 



for - « v 
Ti\av] f v 

for ^«2 <<1 ' 



(B5) 



where the cross section saturates at 



f 2 + (4Am) : 



(B6) 



Accordingly, when the energy scale of the dark matter 
is above the saturation threshold, [<Jv] BV/ (v) oc t>~ 4 , and 
when the cross section is saturated, it is modified by a 
factor of ~ (16v s )~ 2 from the freezeout cross section. 
Behaviors for other cases can be similarly derived starting 
from Eq. (|B2[) . and their features in the extragalactic 
gamma-ray spectrum due to dark matter annihilation can 
then be studied using the methods in this paper. 
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